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By using the concept of the quantum discord (QD), we study the spin-1/2 antiferromagnetic 
Heisenberg chain with next-nearest-neighbor interaction. Due to the SU{2) symmetry and Z2 
symmetry in this system, we obtain the analytical result of the QD and its geometric measure 
(GMQD), which is determined by the two-site correlators. For the 4-site and 6-site cases, the 
connection between GMQD (QD) and the eigenenergies was revealed. From the analytical and 
numerical results, we find GMQD (QD) is an effective tool in detecting both the first-order and 
the infinite-order quantum-phase-transition points for the finite-size systems. Moreover, by using 
the entanglement excitation energy and a universal frustration measure we consider the frustration 
properties of the system and find a nonlinear dependence of the GMQD on the frustration. 

PACS numbers: 03.67.-a, 64.70.Tg, 75.10.Jm 



I. INTRODUCTION 

Quantum correlation is one of the most popular re- 
search topics since it plays a central role in quantum in- 
formation and communication. Usually people thought 
that the quantum computation devices should get their 
computational power from entanglement — one of the 
most essential nonclassical features in quantum mechan- 
ics. Recently, it was found that there exist other nonclas- 
sical relations apart from quantum entanglement. The 
concept of the quantum discord (QD) which is defined 
as the difference between the quantum mutual informa- 
tion and the classical correlation was introduced by 01- 
livier and Zurek [1] to quantify the non-classical corre- 
lations. It has been observed that the QD provided a 
larger region of quantum states with non-classical cor- 
relations, for example, even some separable states have 
non-zero QD [2, 3]. In fact, only zero discord is a nec- 
essary condition for strictly classical correlations [4], so 
that the states with non-zero QD are responsible for the 
efficiency of a quantum computer [5]. Therefore, QD 
could be a new resource for quantum computation, and 
even far cheaper and easier to maintain in the lab. 

Quite recently, many people devoted into the study of 
quantum discord [6-19]. Although there are several ad- 
vantages, a big problem for the application of discord 
is that it is complicated to calculate QD analytically. 
Even for the two-qubit systems, the analytic results of 
QD can only be obtained for a few cases such as X-typc 
states [20], and a general method still lacks. Since that, 
Dakic et al. [21] study the necessary and sufficient condi- 
tion for the existence of non-zero QD for bipartite states 
and introduce a geometric measure of quantum discord 
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(GMQD), which can be evaluated for an arbitrary two- 
qubit state. The concept of the GMQD simplifies the 
calculation of the QD. 

In quantum many-body correlated system, quantum 
phase transition (QPT) is an essential phenomenon and 
attracts widespread attention. As a quantum critical 
phenomenon, QPT happens at absolute zero tempera- 
ture, at which the thermal fluctuations vanish so that 
no classical phase transition could occur. Thus QPT is 
driven only by quantum fluctuation and caused by chang- 
ing the system's Hamiltonian, such as an external mag- 
netic field or the coupling constant. At the quantum 
phase transition point, the QPT behaves as the configu- 
ration transition of the ground state (GS). Therefore, one 
can easily expect that some quantum concepts closely re- 
lated to the ground states can be used to indicate QPTs. 
For example, the concepts of quantum entanglement [22- 
27] , quantum fidelity [28-33] and quantum squeezing [34] , 
have already been widely and successfully employed to 
study QPTs. 

As a concept of quantum correlation, it is natural to 
consider the relation between the QD and the QPTs. Re- 
cently, people become to make use of QD to investigate 
QPTs. In Ref. [35], it shows that QD spotlights the QPT 
point for XXZ Heisenberg chain even at finite tempera- 
ture. On the other hand, the quantum criticalities in the 
environment also play an important role in the dynamics 
of the QD [36, 37]. In some typical systems with QPTs, it 
is found that both classical correlation and quantum dis- 
cord exhibit signatures of the QPTs [38] . All these works 
imply that QD is an effective tool in detecting QPTs. 

In this paper, in terms of the concept of QD, we will 
consider the QPTs of the spin-1/2 Heisenberg chain with 
next-nearest-neighbor interaction which is also called the 
J1-J2 model. This is an interesting quantum many-body 
system for the existence of competition between nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) cou- 
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plings. It well describes the material structure in some 
quasi-one-dimensional compounds, such as CuGeOa [39]. 
There are two important QPTs in this model [40, 41], 
a first-order QPT and an infinite-order QPT. The first- 
order QPT lies on the energy-level crossing of the ground 
state (GS), while the infinite-order QPT is found to be 
closely related to the energy-level crossing of the low- 
lying excited states (ESs) and can be detected by the 
first-ES fidelity [33]. However, the properties of the QD 
in this system is lack of study, which motivates us to 
consider and show that the QD of the GS and the ESs is 
effective to indicate the QPTs in this system. 

This paper is organized as follows: in Sec. II, we in- 
troduce the conception of the QD and the GMQD. In 
Sec. Ill, we give a general analytical result of the GMQD 
for the X-type states. In Sec. IV, we analytically cal- 
culate the eigenenergies of the Ji- J2 model for the 4-site 
and 6-site cases, then obtain the two-site GMQD and QD 
analytically. In Sec. V, we will numerically show that the 
GMQD is a good candidate to indicate the QPTs. In 
Sec. VI, we consider the frustration by making use of the 
entanglement excitation energy (EXE) and a universal 
measure /. Then we find the relation between frustration 
and the GMQD. Finally, the conclusion will be given. 



II. QUANTUM DISCORD AND ITS 
GEOMETRIC MEASURE 



A. Quantum discord 

Given a quantum state p in a composite Hilbert space 
H = Ha "81 "Hb, the total amount of correlation is quan- 
tified by quantum mutual information 



(1) 



where S{p) = — Tr [plog2 p] is the von Neumann entropy 
and pa{B) = T^t^b{A)P is the reduced density matrix by 
tracing out system ]3{A). If we take the system A as the 
apparatus, the quantum discord is defined as follows [1, 
42] 



PB\k = TT^A[iEk<E)T)p]/pk is the postmeasurement state of 
B after obtaining outcome k on A with the probability 
Pk = Tr[{Ek (8) I)p]. {Ek} is a set of one-dimensional 
projectors on Ha, and I is the 2x2 identity operator. 



B. Geometric measure of the quantum discord 

In Ref. [21], Dakic et al. proposed a geometric measure 
of quantum discord defined by 



Dg{p) := mill \ \p- 
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where Q,q denotes the set of zero-discord states and 
WXW^ Tr(XtX) is the Hilbert-Schmidt norm. For 
two-qubit systems, a general state can be written in the 
Bloch representation [21, 43]: 
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i 

with the real parameters Xi, yi, and Rij , and the Pauli 
matrices ai^i^2.3- Here we only consider the case that the 
measurement is taken on the system A. Then an explicit 
expression of the GMQD is obtained as [21]: 







(7) 



where x = {xi,X2,X3)'^ , R is the matrix with elements 
Rij, and femax is the largest eigenvalue of matrix K = 
xx^ + RR^. 

Here we introduce an alternative form which will be 
convenient when we consider the evolution of the GMQD 
under decoherence. First, we introduce a matrix TZ de- 
fined by [44] 



7^ 



1 

X R 



(8) 



and another 3x4 matric TZ' obtained through deleting 
the first row of R: 



D{p)=I{p)-Ca{p), 



(2) 



7^' = (:x,R). 



(9) 



which is the difference of the total amount of correla- 
tion X{p) and the classical correlation Ca{p)- Here the 
classical correlation is defined by 



C(p) = inaxl(pl{ii;fc}), 



(3) 



where X(p]{£'fe}) is a variant of quantum mutual informa- 
tion based on a given measurement basis {Ek\ on system 
A as follows 



Here TZ is just an expectation matrix with the elements 
TZij = Tr[pai (8 aj] for i,j = 0,1,2,3, and ctq = I is 
defined. The expression (9) leads to = 7^'(7^')^. Af- 
ter singular value decomposition, we have TZ' = UAV^, 
where U and V are 3x3 and 4x4 orthogonal matri- 
ces, and A has only diagonal elements A^- = piSij with 
Pi the so-called singular values of the matrix TZ' . Then 
the eigenvalues of the matrix K can be expressed as /J,f . 
Gonsidering [Ja;]]^ -|- [Ji?]]^ = TrK, we get an alternative 
compact form of Dg{p): 



Ap\{Ek}) = S{pb) - ^PkS{pB\k)- 



(4) 



Dgip) 




(10) 
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where the summation and maximization ar( 
all the non-zero singular values /x^ of TV . Thi, 
form will be convenient when we calculate th 
the J I- J 2 model. 



C. The Geometric measure of the quant 
for X-state 
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In the standard basis of operator Sz = S] 
density matrix of the so-called "X-state" is s 
lows: 

(a .9 \ 
_ I b w q\ 

u;* c ' 
\9* d/ 

where the parameters a, &, c and d are real numbers and 
satisfy a+h+c+d — 1, and the positive condition requires 
hc> |u)|^ and ad > \g\^ . Then we can obtain parameter 
matrix R as: 



Rii 


R12 





R21 


R22 











R33 



where Ru = 2Re(g -I- w), R12 = 2Im(w — g), R21 = 
—2lm{g + w), R22 = 2Re('u; — 5), and R33 = a — b — c + d. 
The parameter vectors x = (0, 0, x^)^ with x^ = a + b — 
c—d, and the vector y = (0, 0, ya)^ with j/3 = a — &+c — d. 
By substituting the elements of R and x into the matrix 
7^' = (Eq. (9)), we can calculate the GMQD by 

using Eq. (10) that 

^aiP) = ^ [/^i + Ai2 + Ms - max {^4,^4)] . (13) 

where ^ll = A{\g\ + \w\f, /i^ = 4(|g| - 1^1)' and = 
2 [(a — c)^ -f (5 — d)^] . The expression above is the ana- 
lytical result of the GMQD for the X-statcs, which rep- 
resent an important class of quantum states, e.g., the 
S'C/(2)-symmetric states. 

III. GMQD FOR THE HEISENBERG MODEL 
WITH NNN INTERACTIONS 

The Hamiltonian of one-dimensional Heisenberg sys- 
tem with NNN interaction, i.e., the J1-J2 model reads 
(let h=l) 

N 

H = ^ {JlSi ■ Si+i + J2Si ■ Sj+2) , (14) 

i 

where the denotes the spin-1/2 operator at the i^h 
site. TV is the total number of sites and here we only con- 
sider the even case. Ji and J2 are the nearest-neighbor 
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FIG. 1: Low-level energy spectrum of the spin-1/2 Heisenberg 
chain with NNN interaction for different system sizes = 
4, 6, 8, 10. The parameters are in units of Ji. 



(NN) and next-nearest-neighbor (NNN) exchange cou- 
plings. As usual, we choose the periodic boundary con- 
dition. Note that no exact analytical results are avail- 
able for this model (14) except for the special case of 
J2/J1 = and J2/J1 ~ 1/2. There are two important 
QPTs in this model, a first-order QPT at Jd = 0.5 and 
an infinite-order QPT at Jc2 — 0.241 (where Jc = J2/ Ji)- 

At the critical point Jci = 0.5, the system reduces to 
the Majumdar-Ghosh model [40]. Its GS is of spin zero, 
but degenerate, which is a uniformly weighted superposi- 
tion of the singlet states between NN sites (for even and 
infinite N cases). This point is just the GS energy- level 
crossing induced by the translation symmetry breaking, 
thus it is a first-order QPT point. 

At the other critical point Jc2 — 0.241, the system 
undergoes a Berezinskii-Kosterlitz-Thouless (BKT) type 
QPT from spin fiuid to dimerized phase [41]. This phase 
transition is driven by the competition between the NN 
and the NNN interactions. When J2/J1 < Jc2, the NNN 
interaction does not change the character of the simple 
antiferromagnetic case J2/J1 = 0, whose GS is described 
as spin fluid massless spinon excitations. When J2/J1 > 
Jc2, the frustration term is relevant and the GS flows 
to the strong-coupling dimerized phase. Furthermore, it 
is found that Jc2 is accurately the energy-level crossing 
point of the first ESs for even-size and infinite-size cases. 

As shown in figure 1, wc plot the low-energy levels 
of the GS and the first and second ES versus different 
coupling J2 (let Ji = 1). Clearly, for different sizes of 

= 4, 6, 8, 10, the GS energy-level crossing occurs at 
the point Jd = 0.5, which exactly corresponds to the 
first-order QPT point. Turning to the energy levels of 
the first ES, we find that with the increasing system size, 
the crossing position moves from 0.25 to 0.244, which im- 
plies that in the thermodynamic limit, the energy-level 
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crossing point of the first ES will tend to the BKT QPT 
point Jc2 — 0.241. From the numerical results above, we 
emphasize that the QPTs are closely related to the fun- 
damental change of the energy structure. Consequently, 
the corresponding eigenstates will present special quan- 
tum properties at the QPT points. 

Chen et al [33] studied the ground-state fidelity and 
first-excited-state fidelity of this system. Xiong et al [29] 
found the reduced fidelity can also be used to indicate 
the quantum criticalities. In this paper, we will show 
that the two-site QD of the GS and the first ES is also 
effective to detect the QPTs. 

Due to the SU{2) symmetry and the Z2 symmetry 

generated by commutation [af^^^,H] = 0, the reduced 
density matrix for the two-site subsystem pij =Trjj (ptot) 
with Ttjj meaning tracing out the subsystems except the 
i-th and j-th ones: 



/ a 

6 w 

w b 

VO a 



(15) 



where all the elements are real numbers. Note that if the 
states are degenerate, the reduced density matrices may 
not be in the form above. For example, for the even- 
size J1-J2 Heisenberg system, the first ES is threefold 
degenerate when < J2/J1 < Jc2 {Jc2 — 0.241), and is 
nondegenerate when J2I Ji > Jc2 except for the crossing 
points. In order to overcome the above subtle problem 
induced by the degeneracy, we mix the degenerate states 
with equal probability as that 



Pauli matrices of the two-site subsystem, i.e., (oi^a'^j^p) 
(with a, /3 = x, y, z ). In addition, there is an exchange in- 
variance in the Hamiltonian, which leads to the fact that, 
any NN-site correlators, i.e., (ci^acri+i,^), equal to each 
other, so to the NNN-site correlators {(yi,aO'i+2,p)- Thus 
when we consider the reduced density of the two-site sub- 
system, the elements can be described by the correlators 



(1 + (cr.^CTj^)) , 



(18) 



where cFx.y.z are Pauli matrices. In this system, there 



exists another relation that {(JiaCTja) 
x,y,z and aiaj = a^xO-jx + (^iv<^jy - 



j ) with a = 
Therefore, 



we can prove that the three values fJ.1^2.3 in Eq. (17) equal 
to each other. Then the GMQD becomes 



1 2 



1 

18 



(19) 



This equation is a key result of this paper, and gives a 
general result for the Heisenberg systems which contains 
the same symmetry properties of the J1-J2 model, such 
as XXX-type and dimerized Heisenberg chain. 

Furthermore, for the J1-J2 model, by using of the 
Feynman-Hellman theorem, one can find the correlators 
between the NN sites is 



Ji 



e - J2 



de 



(20) 



(16) 



where G denotes the degeneracy of the energy En and 
\'ipnv) the u-th degenerate eigenstate of the system. 
This assumption is reasonable when wc consider a gen- 
eral mixed state in the thermal equilibrium p{T) — 
exp[— i7/(fcbr)]/Z, with Z is the partition funciton. If 
the eigenstate [-(/'n) of the Hamiltonian H is degenerate, 
each of its degenerate states will have an equal mixture 
weight in the mixed thermal state. It is easy to prove 
that the equal-probability mixed state holds the SU{2) 
and Z2 symmetry, i.e., the reduced density is still in the 
form of Eq. (15). 

By using the results of GMQD for X-statc in Eq. (13), 
we obtain the GMQD for the reduced density matrix pij 
as 



Dg (py ) 



1 r 2 



Ms 



max 



where fif — 4(a — and ^^2=^3= 4?/;^. 

It is known that the elements of the reduced density 
matrix can be presented by the expectation values of the 



where e = (H) /N denotes the expectation value of the 
energy corresponding to each site. Similarly, for the NNN 
sites, we can obtain the correlator {aiai+2) = ^de/dJ2- 
Finally the NN-site GMQD will be connected to the en- 
ergy that 



and the NNN-site GMQD is 

1 



'i,i+2 ) 



.2 8 de 
18 =9 9J2 



(21) 



(22) 



The two equation above present the connection between 
the GMQD and the energy. We can also give the QD of 
this model by using the results of [16] 

D{Pij) = -2(a + 6)log2(a + 6) + 2alog2a + 261og2& 
+2w log2 w-^[{l + 4u7^) log2(l + 4ui^) 
+ (l-4u;2)log2(l-4i«2)], (23) 

by substituting the Eqs. (18) and (20), one can find the 
relation between the QD and the energy. 
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Prom the Eq. (21), we find that the GMQD is domi- 
nated by the energy structure and its derivative. Obvi- 
ously, the GMQD will be sensitive to the singularity of 
the energy level. In other words, the GMQD can effec- 
tively indicate the energy crossings. As is known, the 
first-order QPTs usually lies on the energy-level cross- 
ing of the GS. Hence, we can reasonably expect that 
the GMQD can indicate the first-order QPTs effectively. 
In the J1-J2 model, there is another QPT point nearby 
Jc2 — 0.241, which corresponds to the energy-level cross- 
ing of the first ESs. Therefore, we can use the GMQD of 
the first ES to detect the second QPT point Jc2- 

In the following sections, we will analytically and nu- 
merically calculate the small-size cases, which can also 
justify the fact that the GMQD can indicate the QPTs 
in the J1-J2 model. 



IV. ANALYTICAL RESULTS FOR THE CASES 

OF iV = 4 AND iV = 6 

In this section, we will analytically calculate the Ji- 
J2 model for 4-site and 6-site cases. The eigenenergy 
structure can be directly used to obtain the GMQD. In 
the following, we set Ji = 1 for simplicity. 

Due to the periodic boundary condition, the Hamil- 
tonian is translational invariant, i.e., [H^T] = 0, 
where T defines the cyclic right-shift operator satisfying 
T |mim2...TOi-i7Tii) = \mimi...mi-i) . 



A. 4-site case 

For [H, Jz] = 0, the 16-dimcnsional Hilbert space for 
a 4-site J1-J2 model can be divided into invariant sub- 
spaces spanned by vectors with a fixed number of re- 
versed spins [45]. Thus, the dimension of the largest sub- 
space is 6. With the help of translational invariance, we 
can further reduce the Hamiltonian matrix to 2 x 2 sub- 
matrices, and then the eigenvalue of Hamiltonian can be 
solved. 

The none-reversed subspace contains only one vec- 
tor 1 0000), which the eigenvalue is i? = 1. The sub- 
space with one reversed is spanned by 4 basis vectors, 
{T" 1 1000) , n = 0, 1, 2, 3}. Considering the translational 
invariance of the Hamiltonian, we choose another set of 
basis vectors [46] : 

3 

=5^c.^T"|1000), (24) 

where ujk = e''^^'''^,k € {0, 1,2,3}. One can check that 
I'^k) is eigenstates of T with eigenvalues uj^^, and also are 
eigenstates for H with eigenvalues i (— 1 + ^ + J2i^k^ + 

When the reversed number is 2 or 3, we choose the 



basis 

3 

|vl/fe) =^c.,"r" 11100), (25) 

n=0 

and 

|<I>,)=^c^^r" 11010), (26) 

n=0 

for the 6-dimensional subspace. Under this space, one 
can rewrite the Hamiltonian into a 2 x 2 form and obtain 
all the eigenvalues. Then, the GS energy is 

( -2 + J2 when J2 < 0.5, , . 

\ -3J2 when ,/2> 0.5. ^'^'^ 

Substituting the energy levels above into Eq. (21), one 

can obtain the NN-site GMQD corresponding to the GS 
Po- 

We also obtain the 1st ES energy: 

{-I + J2 when J2 < 0.25, 
-3J2 when 0.5 > J2 > 0.25, (29) 
-2 + J2 when J2 > 0.5, 

thus the NN-site GMQD corresponding to the first ES 
Pist: 

( ^ when J2 < 0.25, 
Dgipist) = < when 0.5 > J2 > 0.25, (30) 
[ I when J2 > 0.5, 

Obviously, the GMQD presents discontinuity at the point 
J2 = 0.5, which just corresponds to the first-order QPT 
of this system. While for the first ES, the GMQD sudden 
changes at two points: J2 = 0.25 and J2 = 0.5. Similarly, 
one can obtain the NNN-site GMQD, and we do not show 
the results here. 



B. 6-site case 

For the 6-site case, the energy structures are more 
complicated. Fortunately, we analytically obtain all the 
eigenenergies. The GS energy is: 

^ f -5 V9J| - I8J2 + 13-1 when J2< 0.5, 
1 -f(l + J2) when J2> 0.5. 

(31) 

Therefore, when J2 < 0.5 the NN-site GMQD for the GS 
is 
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(a) (b) (a) (b) 



FIG. 2: (a) NN-site GMQD of the ground state versus cou- 
pling J2 for different system sizes A'^ = 6, 8, 10; (b) The 
derivative of the GMQD versus J2. The parameters are in 
units of Ji. 



FIG. 3: (a) NN-site GMQD of the first excited states versus 
coupling J2 for different system sizes A'^ = 6, 8, 10; (b)The 
derivative of the GMQD versus J2. The parameters are in 
units of Ji. 



I8J2 - 26 



162 \ V9J| - I8J2 + 13 



(32) 



and when J2 > 0.5, we have a constant value of the 
GMQD: Dg{pQ) = 1/18. 
The first ES energy is: 



and the critical values is obtained in the systems with 
infinite size. Although the small systems cannot accu- 
rately reflect the quantum criticalities, we can simulate 
the critical properties in the small systems [46, 47]. In 
the following section, we will numerically calculate the 
GMQD for some larger system with 8 sites and 10 sites. 
The second critical region can be approximately indicate 
by the GMQD of the ES. 



-i^gJl- IOJ2 + 5- 1 when J2 < 0.25, 



E = { -|(1 + J2) when 0.5 > J2 > 0.25, 

-iv^Q-^l -I8J2 + I3-I when J2> 0.5. 

(33) 

Thus when J2 < 0.25, the NN-site GMQD is 




9J2 



19 J2 + 10 



2 V9J| - IOJ2 + 5 



(34) 



When 0.5 > J2 > 0.25, we have a constant value of 
the GMQD that Dg{pist) = 1/18. When J2 > 0.5, the 
GMQD becomes 



Dgipist) = 



9J| 



81 



27J2 + 22 



2v/9J| - 18 J2 + 13 



+ 1 



(35) 



Obviously, the GMQD of the GS changes at the point 
J2 = 0.5, and the GMQD of the ES changes at both the 
points J2 = 0.5 and J2 = 0.25. 

From the analytical results of the GMQD, we see that 
the first-order QPT point Jci =0.5 can be indicated by 
the GMQD of the GS and ES. However, for the small- 
size cases as 4 sites and 6 sites, the GMQD of the ES 
cannot indicate the second QPT point Jc2 — 0.241. This 
is because the QPTs occurs at the thermodynamic limit 



V. NUMERICAL RESULT 

In Fig. 2a, we plot the NN-site GMQD of the GS versus 
coupling J2 for different system sizes iV = 6, 8, 10. It is 
seen that the GMQD presents discontinuity at the point 
Jci = 0.5. This is because that the structure of the 
GS changes suddenly at the QPT point which can be 
refiected by the quantum correlation, i.e., the GMQD. In 
Fig. 2b, the QPT point Jd = 0.5 is clearly detected by 
the derivative of the GMQD. 

In order to study the BKT-type QPT, in Fig. 3a we 
plot the NN-site GMQD of the first ES versus coupling J2 
for different system sizes N = 6,8, 10. Obviously, there 
exists a sudden drop of the GMQD. Moreover, with the 
increasing system size, the point of the drop tends to the 
QPT point Jc2 — 0.241. In Fig. 3b, it is more obviously 
to find the QPT point when we plot the derivative of the 
GMQD. The NNN-site GMQD of the first ES is shown 
in Fig. 4. Some different phenomena can be seen, for 
the case A^ = 6, increasing coupling J2 can enhance the 
NNN-site GMQD, and the derivative of the NNN-site 
GMQD presents a positive peak. With the increasing 
system size, the discontinuous position of the NNN-site 
GMQD approaches the QPT point Jc2 ^ 0.241. 

Note that as the system size increases, the discontin- 
uous behavior around the QPT point Jci = 0.5 and 
Jc2 — 0.241 becomes more and more weak. One pos- 
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0.24 0.245 0.25 0.255 0.26 
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2 0.22 0.24 0.26 0.28 0.3 
J, 



FIG. 4: (a) NNN-site GMQD of the first excited states versus 
coupling J2 for different system sizes N = 6, 8, 10; (b) The 
derivative of the GMQD versus J2. The parameters are in 
units of Ji. 



siblc explanation is that for the QPTs associated with 
the energy-level crossings, the continuity of the quantum 
correlation characterized by the GMQD lies on the or- 
thogonality of the eigenstates on both sides of the point. 
However, in this paper we only consider the GMQD in 
the reduced two-site subsystem where the orthogonal- 
ity of the global system is destroyed. Moreover, with 
the increasing system size, the orthogonality, in other 
words, the information kept by the two-site subsystem 
occupies less and less proportion of the global system. 
Therefore, the discontinuous phenomena of the GMQD 
become weaker and weaker with the increasing system 
size. Nonetheless, this GMQD approach for small finite- 
size systems is still meaningful for the usual theoretical 
and experimental researches in the QPT problems. 

We also numerically calculate the QD (defined in 
Eq. (2)) and find a similar discontinuous behavior with 
the GMQD near the QPT points, thus we do not show 
the numerical results of the QD in this paper. It is known 
that the GMQD can not absolutely characterize the QD, 
e.g., under a local decoherence channel, a sudden change 
in the decay rate of the GMQD does not always imply 
that of the quantum discord [44]. However in this J1-J2 
model, the GMQD reflects the quantum criticalities as ef- 
fectively as the QD. Furthermore, the concise expression 
of the GMQD is more conducive to reveal the connection 
between the quantum correlation and the eigenenergies. 
The numerical calculation of the GMQD is more time- 
saving than that of the QD. 



VI. FRUSTRATION, GMQD AND QPT 

As is konwn that the J1-J2 model characterizes a typi- 
cal frustrated spin system. Frustration arises from the si- 
multaneous presence of competing antiferromagnetic ex- 
change interactions of different spatial range [48] , and 



which causes the impossibility of minimizing simultane- 
ously the energy of competing interactions. Some meth- 
ods based on quantum information techniques were in- 
troduced to investigate frustrations, e.g., the entangle- 
ment excitation energy (EXE) AE [49] , and a universal 
measure of frustration / [50]. In this section we will an- 
alytically calculate the two quantities and find the rela- 
tionship with GMQD. 

The entanglement excitation energy (EXE) is defined 

as 



= min ((*[ UkHUk [*) - (*| H [*)) , 



(36) 



where Uk = ®i^kli ® Ok with li denoting the identity 
operator on all the spins but the one at site fc, and Ok is 
a generic Hermitian, unitary, and traceless operator can 
be described as Ok = smO cos ipak,x + sinO sunpok^y + 
cosOcTk.z- If choosing the as the GS the system 
H, we know that for any translationally invariant and 
frustration-free Hamiltonian H such that [i/, Uk] 7^ 
VC/fe, the vanishing of the EXE is a necessary and suf- 
ficient condition for GS factorization [49] , and thus it is 
a proper measure of single-site entanglement. Indeed, 
the presence of frustration tends to enhance correlations 
among the constituents and thus to depress the possibil- 
ity for the occurrence of separable (uncorrelated) states. 
Hence entanglement and separability can be used to qual- 
ify and quantify frustration [48] . 

When calculating the EXE in the J1-J2 model, due to 
the SU (2) and Z2 symmetry and also the translationally 
invariant, we have the relations {(Ji^aO'j^p) = for a ^ (3 
andQ!,/3 = x,y,z, and {a^^^aj^^) = {(Ji,y(yj,y) = {(Ji,z(y],z) ■ 
Finally we obtain that 



AS,, 



-2Ji ((T; zfJi+i,, 



2J2 



'i,z<Ji+2,z/ 



(37) 

where = (H) /N denotes the expectation value of the 
energy corresponding to each site. Especially in this 
model, no matter what values of the parameters 9 and 
ip are chosen, the EXE AEi gives the same value. Obvi- 
ously, AEi can be used to characterize the QPTs in the 
J1-J2 model due to the direct relation with the energy 
level. The energy level crossing point also corresponds 
to the sudden change of the EXE and thus indicate the 
QPTs. The relation in Eq. (37) is still true in the excited 
states, where AEi may present negative values. There- 
fore, AEi is effective to detect the second critical point 
J2C ~ 0.241 when taking into account the first excited 
state energy. Recalling the GMQD of the NN spin pair 
in Eq. (21), then we have 



and the GMQD of the NNN spin pair 



(38) 



Dgip 



13) 



1 aA^ 2 
8^ dJ2 ' ' 



(39) 
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In the J1-J2 model, the GS of hs = Jis^ • s^+i (or 
J2Si • Si+2) is the singlet state |*_) ^ (|10) - |01)) /V2, 
thus the projector Us = |*-> (*-| with d = 1. Then 
we obtain the frustration measure is a frmction of the 
correlator as: 



Jij ^ ' ^ \'^ia<^ ja) , 



(42) 



where a = x,y^ z, and j = i + 1 or i + 2. 

From the frustration measure in Eq. (42), we exploit a 
nonlinear relation with the GMQD that 



(43) 



FIG. 5: (a) NN-site frustration measure fi.i+i and its lower 



bound eI^^^-^ versus coupling J2 for the system size TV — 10; 
(b) NNN-site frustration measure fi.i+2 and its lower bound 
versus J2 for the system size A'^ = 10. The coupling 
strength J2 is in units of Ji. 



the two equations above connect the pairwise GMQD to 
the EXE and its derivative. A vanishing EXE will induce 
a zero pairwise GMQD but the converse may be not true. 
Finite EXE can also cause vanishing GMQD, because 
GMQD characterizes the pairwise quantum correlations 
and not the single-site ones. In other words, frustration 
can destroy the pairwise correlation and hold the block 
correlation between single-site and the rest part. 



A. Universal measure of frustration and GMQD 

Recently a universal measure of frustration has been 
introduced as the overlap between the global ground state 
of a frustrated model and the GS of unfrustrated spin 
pairs [50]. For a many-body Hamiltonian H = J^s^s 
with hs describing the local interaction. The definition 
of the frustration corresponding to the subsystem S is 



fs = l- TApUs 



(40) 



where p denotes a pure state, and let Us = |'I')^ (^"1 be 
the projector onto the subspace S and 1^')^ is the GS 
of the local interaction hs which is not frustrated. The 
quantity fs is a well-defined measure of the interaction 
hs- It turns out that the geometric bipartite ground- 
state entanglement measured by 



E 



id) 



1 - 



At {Ps 



(41) 



is a universal lower bound to the frustration fs, and d 
is the rank of the projector Us- For the degenerate GS 
d indicates the degeneracy. are the eigenvalues 
of Ps =Trij (p) in decreasing order. 



from which, the GMQD achieves the maximum value 1/2 
when the frustration vanishes fij = 0, and has zero value 
when a finite frustration exists fij ~ 3/4. This implies 
that the frustration in this model will depress the pair- 
wise correlation, however, it enhances the block correla- 
tion between the two-site subsystem and the rest sub- 
system, which can be verified by calculating the linear 
entropy 



5Kp.,) = l-Tr(p2 



---- f -- 
- 4 3 V'' 4 



(44) 



obviously, when fij = 3/4, the linear entropy has the 
maximum value which signals the the maximal bipartite 
entanglement between the two-site subsystem and the 
rest part of the system. To note that we now consider the 
pure GS case, thus the linear entropy is an effective mea- 
sure of the bipartite entanglement and also the quantum 
correlation because we also find a direct relation between 
the linear entropy and the GMQD that 



Si {Plj) = ^ - l^Dg (Pij) 



(45) 



from which one can clearly find that the maximal pair- 
wise correlation Dg [pij) = 1/2 will destroy the block 
correlation between the two-site subsystem and the rest 
part, then Si (pij) = 0. 

The total measure of frustration is defined as _F = 
XT J fij J where M is the number of bonds, and in our 
model M — 2N (N is the total number of the spin), then 
we obtain 



F 



2 + fi:i+'2) 

4^8 ^ {'^i.a<^i+2.a)) 



3 1 



[Ji - J2 



in the above, we make use of iuiniyi 



^2!^)' and 



4 

3,7i 



(46) 



(e. 



.en I 



4 

3 aj2 



corresponding to the 



NN sites and NNN sites respectively. 
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Now we calculate the lower bound 



4' - u 



when 



<o, 



Ciq) when ((TiQCTj-Q, ) > 0. 



(47) 



With the help of the energy levels shown in Fig. 1, for the 
GS energy level, obviously, one can find (criQ:Cri+2.a) > 
when Ji < 0.5, and {aiaa'i+2,a) < when J2 > 0.5. Thus 
based on the Eq. (47), we have 



fi,t+2 > £^,v,V2 '^hen J2 < 0.5, 



■^hen J2 > 0.5. 



(48) 



It is known that the inequality fij > Efj implies the 
appearance of the geometric frustration [50] . Hence, we 
find that the geometric frustration appears when J2 < 0.5 
and vanishes when when J2 > 0.5. It implies that the 
frustration measure and its lower bound can reflect the 
change of the fundamental structure of the system, such 
as the energy level crossing, therefore, they can be used 
to indicate the first-order QPT point Jd = 0.5. 

Based on the analytical results of the GS energy level 
in the 4-site and 6-site case, we calculate the frustration 
measure. For 4-site case (let Ji = 1): the frustration of 
the NN sites 



fi,i+l — 



1 = when J2 < 0.5, 

I = ^M+i when J2 > 0.5. 



(49) 



and the frustration of the NNN sites 




= £;|V+2 when J2 > 0.5. 
For 6-site case, we have the frustration of the NN sites 
= E^, when J2 < 0.5, 



^ 1 J= E^l, when J2 > 0.5, 



(51) 



where O = ■\/9J| — I8J2 + 13, and the frustration of the 
NNN sites 



3(0-./2 + l) . p(i) _ 30+./2-1 

1 - ^(1) 

2 ~ ^i,i+2 



_j 4^ > 41^2 = ^^^^^ When J2 < 0.5, 



\ = e\^},^ when J2 > 0.5. 



(52) 



For larger size system such as = 10, we numerically 
calculate the frustration measure and its lower bound 
in Fig. 5. Both the NN-site frustration measure 
in subfigure (a) and the NNN-site frustration measure 
fi.i+2 in subfigure (b) present a sudden change at the 
QPT point J2 = 0.5. Moreover, in the subfigure (a), we 
shows that the NN-site frustration measure fis+i and its 
lower bound E^^^j^-^ are consistent with each other, which 
implies that only quantum frustrations exist between NN 



sites. Differently, in the subfigure (b), we find the NNN- 
site frustration fi^i+2 > £'lV+2 when J2 < 0.5, which 
signals the appearance of the geometric frustration. 

From the above, we find the nonlinear dependence 
of the GMQD on the frustration, which implies that 
the pairwise quantum correlation characterized by the 
GMQD is greatly affected by the frustration, e.g., some 
finite frustration as =3/4 can depress the GMQD 
to zero. On the other hand, we believe that the frustra- 
tion measure can be used to characterize the first-order 
QPT in the J1-J2 model and also other systems. In ad- 
dition, the inequality fs > Eg'''' holds as well in any 
mixed states [50] , thus the method is also effective to de- 
tect the first-order QPTs in the degenerate GSs. In our 
J1-J2 model, the quantity fs can also be used in the 
first-excited states to detect the QPT point Jc2 ~ 0.241. 
However, in the excited states (ESs) case, the quantity 
fs can not be understood as the frustration measure, in- 
stead, it only quantifies how much fails to fully overlap 
with the subspace selected by the projector Hg, and the 
projector H5 may come from the GS or the ES of the 
local interaction ^15. 



VII. CONCLUSION 

In this paper, we considered the quantum discord (QD) 
in the Heisenberg spin chain with next-nearest-neighbor 
(NNN) interaction. By using the geometric measure of 
the quantum discord (GMQD), we studied the quantum 
correlation properties of the ground states (GSs) and the 
first excited states (ESs). 

We give a general analytical result of the GMQD for 
the X-type states. For the Heisenberg system with the 
SU (2) symmetry and Z2 symmetry, we give an exact 
relation between the GMQD and the two-site correla- 
tors. Furthermore, the connection between GMQD and 
the eigenenergies was revealed. For the 4-site and 6-site 
cases, the analytical results of the GMQD for the GS and 
first ES are obtained, from which the first critical point 
Jci = 0.5 can be exactly detected. We also numerically 
studied the NN-site and NNN-site GMQD of the ESs. It 
is found that when the system size increases from 6 sites 
to 10 sites, the discontinuous point of the GMQD tends 
to the second QPT point Jc2 ^ 0.241. 

Moreover, by using the entanglement excitation energy 
and a universal frustration measure we considered the 
frustration properties of the system and find the non- 
linear dependence of the GMQD on the frustration. The 
measure of the frustration can also be employed to detect 
the QPTs in this system. 

We emphasize that the two-site GMQD can detect the 
QPTs in this Heisenberg system with NNN interaction. 
Although the two-site GMQD approach is only effective 
in the finite-size systems, it has practical significance 
for the usual theoretical and experimental studies. The 
problems of the GMQD (or QD) in other systems with 
QPTs are interesting and need further consideration. 
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